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Abstract 

A continuum model of dislocation pileups that takes the self-energy of dislo¬ 
cations into account is proposed. An analytical solution describing the distri¬ 
bution of dislocations in equilibrium is found from the energy minimization. 
Based on this solution we show (i) the existence of a critical threshold stress 
for the equilibrium of dislocations within a double pileup, and (ii) the exis¬ 
tence of a non-linear regime in which the number of dislocations in a double 
pileup does not scale linearly with the resolved external shear stress, contrary 
to the classical double pileup model. 
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1. Introduction 

In recent years there is a substantial amount of literature dealing with the 
dislocation pileups in crystals within the continuum approach (see, e.g., mu 
[3 DSl DSl El UHl DSl [201 EH 1221 EH and the references therein). The proposed 
models turn out to be quite useful as they can predict the size effect for the 
yield stress, which agrees quite well with the experimental data 121 but 
however does not conhrm the well-known empirical law formulated by Hall 
[9] and Fetch [25] . A traditional explanation of the Hall-Fetch relation, based 
on the classical dislocation pileup model considered in [20] , is that dislocation 
pile-ups serve to enhance the stress felt at grain boundaries. However, the 
Leibried’s dislocation pileup model differs from the contemporary continuum 
models [21IIllI21ll3lEllIHlll9lEniEIlE2] in two important aspects: hrstly 
the absence of a threshold stress for dislocation nucleation, and secondly the 
absence of a hnite-sized dislocation-free region. This is the motivation for us 
to reconsider the Leibfried’s model in order to resolve this discrepancies. 
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Leibfried’s one-dimensional theory of dislocation pileups is based on the 
force equilibrium: the Peach-Koehler resultant force acting on a dislocation 
produced by applied external stress and by other dislocations must vanish 
|3l [23] . Within the continuum approximation one can deduce from here the 
well-known integral equation |23], provided the dislocation density is not 
zero. This equation need not be satished in a dislocation-free zone. In his 
now classical paper Leibfried [23] mentioned that the natural way to hnd 
the stable equilibrium distribution of dislocations if such a zone occurs is 
to use the variational principle of minimum energy: among all admissible 
distributions of dislocations the true stable distribution minimizes energy 
of the crystal (cf. with the LEDS-hypothesis formulated in [10] )• To the 
best of the author’s knowledge, the one-dimensional continuum model of 
dislocation pileups based on the energy minimization taking the self-energy 
of dislocations into account has not yet been proposed. The aim of this 
short paper is to £11 this gap. We start from the expression for the energy 
of crystal containing an array of dislocations which includes also the self¬ 
energy of dislocations. The self-energy density is proportional to the absolute 
value of the dislocation density, so this does not change the resultant force 
acting on the dislocation as well as the integral equation except the forces 
acting at the tails of the pileups. However, this self-energy influences the 
stable equilibrium distribution of dislocations essentially. We will show that 
the density of dislocations is identically zero if the applied stress is found 
below some critical value called the yield (or threshold) stress. Besides, 
the number of dislocations depends on the applied stress non-linearly, in 
contrast to the classical theory. The results of the proposed theory agree 
quite well with those of the continuum models of dislocation pileups obtained 
in [H 1121 HSl HU EHl Ell E2] • A continuum model of dislocation pile-ups taking 
into account the Frank-Read source proposed in [1] (see also H) leads to the 
similar results although it does not have the energetic structure. 

In the next Section we present the continuum model of dislocation pileup 
accounting for the self-energy of dislocations. In Section 3 the energy min¬ 
imization problem is solved and the related property of the solution is dis¬ 
cussed. Finally, Section 4 concludes the paper. 

2. Continuum model of dislocation pileup 

Consider the plane strain problem of an infinite crystal which is uniformly 
loaded by a shear stress r applied at infinity. Under this loading condition a 
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Figure 1: Double pileup of edge dislocations 


linear array of equal number of positive and negative edge dislocations may 
occur on the slip line which is chosen to be the x-axis (see Fig. [^. The 
dislocation lines are parallel to the ^-axis, while their Burgers’ vectors are 
directed along the x-axis. We assume that there are two obstacles (like two 
inclusions or grain boundaries) at x = ±c so that dislocations are confined to 
stay in the interval L = (—c, c) of the x-axis. In the continuum limit we may 
replace the sum of many closure failures induced by dislocations in form of 
step functions by a smooth function (p(x) (see Fig. |^. The obstacles impose 
the following constraints on this function 


(p(±c) = 0. 


( 1 ) 


Now we present the resultant inverse plastic distortion im in the form 


-f^xy = ^{x)5{y), 


with 5{y) being the Dirac-delta function. Differentiating this equation with 
respect to x we obtain 

- I3xv,x = ^'{x)5{y). (2) 

The interpretation of ([^ is quite simple: if we integrate this equation over a 
circle C with the middle point at x = —c and the radius c + x, then 


dxdy 


'c 



^'{0 di = (p(x) 


Thus, we get the closure failure of an amount (p{x) which should be equal 
to the net Burgers’ vector of all dislocations within the interval (—c,x). 
Denoting the dislocation density by p(x), we get 

= b j p(0 d^ => (p'(x) = bp{x). (3) 
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So, it is extremely useful to think of this double pileup of dislocations as a 
mode II cracl^and to interpret the closure failure V3(x) as the crack opening 
which is shown schematically in Fig. 


(p(x) 



Figure 2: Continuum approximation of the closure failure 


The energy of this system can be found directly from the energy of crystal 
containing continuously distributed dislocations in the interval (—c, c) of the 
x-axis. Indeed, for the plane strain state the energy functional per unit depth 
of the crystal reads 

/[u(a;,t/)]= / (l){eap - €li^)dxdy - / TaUads, 

J A JdA 

where A is the cross-section area of the crystal in form of a cylinder, 
the free energy density, the symmetric part of the plastic distortion, 
while Ta is the external traction acting at the boundary dA. Substituting 
0 = ^o'a/ 3 (Ma ,/3 — Pajs) luto this energy functional and integrating the term 
by parts using the equilibrium equation and the boundary conditions, 
we can show that this term is canceled out with the last term, so 




^xyPxy dxdy 


a. 


xy 


[x)ip{x) dx. 


Let us hrst exclude the self-stress (causing the self-energy of dislocations) and 
neglect the influence of the boundary of crystal by assuming that A occupies 
the whole (a:, ?/)-plane. Then we can present the shear stress a^yix) on the 


^In fact, various crack problems have been solved within this continuum model (see, 

for example, mm)- 
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a:-axis in the form 


(^xy\X, 




i-x 


di + r, 


D = 




27r(l 


(4) 


The hrst term in this formula gives the stress held induced by the continuous 
distribution of dislocations excluding those in the neighborhood of x (the 
integral ^ in (|^ is dehned as Cauchy’s principal value), the second term 
corresponds to the stress held caused by the external shear stress r. Substi¬ 
tuting (|^ into the energy functional and integrating the hrst term by parts, 
we obtain hnally 





In \x — ^\ip'{x)(p'{^) dx — / T(p{x)dx. 


(5) 


L JL 


The double integral in ([^ is the symmetric form of the quadratic functional 
which should be understood in the following sense 



ln|a;-^|(p'(a;)<^'(0 dxd^ = 


L JL 


L£ 


y'(0 

^-x 


d^ip{x)dx. (6) 


It is interesting to note that this quadratic functional is positive dehnite and 
obeys an inequality similar to that of Wirtinger (see [IS]). Note also that, 
since the left-hand side is invariant with respect to the scaling x ^ x/c due to 
the boundary condition ([^, the arguments in the log-kernel can be regarded 
as dimensionless. 

Functional ([^ is the energy functional of the classical continuum dislo¬ 
cation theory in which, as we have seen from the previous derivation, the 
self-stress of dislocations is excluded. This means that the self-energy of dis¬ 
locations is neglected. Let us now take the self-energy of dislocations into 
account by the following deliberations. It is well-known that the self-energy 
of one dislocation (per unit length) does not depend on its sign and equals 
K/i6^, where k = ln(4i?/6)/(47r(l — u)), with R being the distance of dislo¬ 
cation to the boundary of crystal and u the Poisson ratio m- Since the 
dislocation density is p{x) = j;(p'(x) and since the self-energy does not de¬ 
pend on the sign of dislocations, the total self-energy of dislocations is given 
by 

Kpb / \p'{x)\dx. 
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Adding this energy to the functional ([^ we get dually 


= -f 



In \x — ^\ip'{x)ip'{^)dxd^— / Tip{x)dx + K^b / \i~p'{x)\dx, 


L JL 


(7) 

We formulate the following variational principle: the stable equilibrium dis¬ 
tribution of dislocations corresponds to the minimizer (p{x) of functional Q 
under the constraints Q. The existence of minimizer is guaranteed by the 
above mentioned inequality for the quadratic functional ([^ provided in [I6] . 


3. Energy minimization 

Symmetry of the problem implies that (p{x) is even and (p'{x) is odd. It 
can be shown that the presence of \(p'{x) \ in the last term of the energy func¬ 
tional Q causes the minimizer to have the dislocation-free zone as observed 
in [1]. Therefore we assume that there is a length a < c such that 

{ >0 for X G T_ = (—c, —a), 

<0 for X G T+ = (a, c), (8) 

0 otherwise. 


Thus, there are two boundary layers where positive and negative dislocations 
pile up against the obstacles and the dislocation-free zone in-between. Let 
us drst £x the constant value of (p{x) in (—a, a), (p(a), as well as the length 
a. Varying functional Q and taking into account Q we get for the drst 
variation 


SI 



^-X 


di-D-l 
J L+ 


^ — X 


d^-T 


6ip{x) dx. 


There is no contribution of the last integral in ([^ due to the fact that it is 
equal to 2Kfib(p{a) which is dx. Since 6(p{x) may be chosen arbitrarily for 
X G ip{x) should satisfy the singular integral equation 




for X G L±. 


(9) 


Now, it is not difdcult to show that the variation of a leads to the boundary 
conditions 


ip\±a) = 0 , 


( 10 ) 
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which guarantee the continuity of the dislocation density. Thus, equation 
is subjected to the boundary conditions ([^ and (10). The solution of (|^, 
([^, and ( 10 ), found in |23] (see also [261 E]), reads 


ip'{x) = 



for a: e L±, 
otherwise. 


From here one easily finds (p{x) in terms of the elliptic integrals (see [5], p 
276) 

^[cE{X,q) - ^F{X,q)] for a; G L±, 

^[cE{q) - ^K{q)] otherwise. 


V7(x) = 


( 11 ) 


where 


A = arcsin 


— x‘^ 

— a^‘ 


q = 


\/c2 - 


In formula ( [IT| E{X,q) and E(X,q) are the incomplete elliptic integrals of 
the first and second kind, respectively, while K{q) and E{q) are the complete 
elliptic integrals of the first and second kind. Mention that q){a)/h gives the 
total number of dislocations of either sign in the pileups (cf. (|^). 

Up to now (^(a) is an unknown quantity. It can be found by minimiz¬ 
ing energy as function of a. Substituting the solution (0 into the energy 
functional ([^ and using ([^ we obtain 


r 


/(r,a) = -- 


(p(a:) dx + 2fiKbip{a) = 


-T / ip{x) dx — Tq:>{a)a +2^Kbip{a). 
Jl+ 

( 12 ) 

The direct consequence of (12) is that for r < 2^Kb/c the minimum of energy 
is achieved at the end-point a = c giving 1 = 0. Indeed, for a < c we have 
(p{x) < (p{a), so, I{T,a) > 0 in this case and the minimum is achieved at 
a = c. Thus, the threshold stress Tc = 2qiKb/c. Note that this threshold 
stress is inversely proportional to the size of the grain exhibiting clearly the 
size effect (cf. with 


In general the integral in (12) cannot be found in closed analytical form. 


However it is easy to find minimum of (12) with respect to a numerically. 


For this purpose let us introduce the following dimensionless variable and 
quantities 

X —I 

C = ®(C) = V?(Cc)/c, a = a/c, (3 = b/c, 7 = r//i, I = —. 

C fjLC^ 
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Figure 3: Graph of function a) bold line: theory with self-energy of dislocations, 

b) dashed line: classical theory 


Then function ‘h(C) is given by 


^(C) = 27(1 - J^)[E{X, q) - a^F{X, q)] for C e (a, 1), 


where 


A = arcsin a / ^q 


l-a2’ 


= a/1 — 


Formula (12) becomes 


1(7,0;) = —7 / <F(C) (iC ~ + 2 fi:/ 3 <F(Q;). 


We evaluate this function numerically using Mathematica. The minimum 
of with respect to a is also sought numerically. After finding a* at 

which energy reaches the minimum, we find the total number of dislocations 
of either sign in the pileups by (cf. ([^) 

N = ^ _ ayr(,(a.))|. 

In Fig. 1^ the graph of N versus the shear strain 7 is plotted, where, for 
comparison, the dashed straight line N = 27(1 — v)/(5 obtained by the clas¬ 
sical theory is also shown. For the numerical calculation we took (3 = 0.001, 
K = 1.0, V = 0.25. It is seen that for 7 < 2Kf3 no dislocations can be 
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Figure 4: Normalized dislocation density p{Q: a) bold line: theory with self-energy of 
dislocations, b) dashed line: classical theory 


formed. The curve A^( 7 ) obtained from the energy minimization is a non¬ 
linear function of 7 , but for 7 > Gk/S it approaches quickly the straight line 

N = 27(1 - u)/f3. 

Fig.0 shows the normalized dislocation density p = nDcp' /t as function 
of the dimensionless coordinate C, = x/c for 7 = 0.003 which, in contrast 
to the classical distribution —(/yjl — (presented by the dashed line), de¬ 
pends on the stress (or the shear strain 7 ) through the half-length of the 
dislocation-free zone a. However, this haft-length approaches quickly zero 
for 7 > O/t/d, so the distribution of dislocations approaches quickly that of 
the classical theory at large 7 . 

4. Conclusion 

It is shown in this paper that the account of self-energy of dislocations 
leads to the existence of the threshold stress at which the dislocations nucle¬ 
ate. The total number of dislocations depends non-linearly on the applied 
shear stress, however it approaches quickly that of the classical theory for 
the stress three times larger than the threshold value. 
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